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Abstract 



Fermionic condensate and the vacuum expectation values of the energy-momentum ten- 
sor are investigated for a massive spinor fields in higher-dimensional spacetimes with an 
^ I arbitrary number of toroidally compactified spatial dimensions. By using the Abel-Plana 

summation formula and the zeta function technique we present the vacuum expectation val- 
ues in two different forms. Applications of the general formulae to cylindrical and toroidal 
^ I carbon nanotubes are given. We show that the topological Casimir energy is positive for 

\^ I metallic cylindrical nanotubes and is negative for semiconducting ones. The toroidal com- 

04 ■ pactification of a cylindrical nanotube along its axis increases the Casimir energy for metallic- 

' type (periodic) boundary conditions along its axis and decreases the Casimir energy for the 

semiconducting- type compactifications . 

(N 
O 

a^ 
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1 Introduction 

Many of high energy theories of fundamental physics, including supergravity and superstring 
theories, are formulated in spacetimes having compact spatial dimensions. From an inflationary 
point of view universes with compact dimensions, under certain conditions, should be consid- 
ered a rule rather than an exception [1]. The models of a compact universe with non-trivial 
topology may play an important role by providing proper initial conditions for inflation. There 
are many reasons to expect that in string theory the most natural topology for the universe is 
that of a flat compact three-manifold [2] . The quantum creation of the universe having toroidal 
spatial topology is discussed in [3] and in references [1] within the framework of various su- 
pergravity theories. An interesting application of the quantum field theoretical models with 
non-trivial topology of spatial dimensions recently appeared in nanophysics [5j. In a sheet of 
hexagons from the graphite structure, known as graphene, the long-wavelength description of 
the electronic states can be formulated in terms of the Dirac-like theory of massless spinors in 
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3-dimensional spacetime with the Fermi velocity playing the role of speed of light (see, e.g., 
Refs. p]). Single- walled carbon nanotubes are generated by rolling up a graphene sheet to form 
a cylinder and the background spacetime for the corresponding Dirac-like theory has topology 
R'^ X S^. Compactifying the direction along the cylinder axis we obtain another class of graphene 
made structures called toroidal carbon nanotubes with the background topology x (5^)^. 

The compactification of spatial dimensions leads to a number of interesting quantum field 
theoretical effects which include instabilities in interacting field theories [7], topological mass 
generation [8], symmetry breaking [9J. In the case of nontrivial topology the boundary condi- 
tions imposed on fields give rise to the modification of the spectrum for vacuum fluctuations 
and, as a result, to the Casimir-type contributions in the vacuum expectation values of phys- 
ical observables (for the topological Casimir effect and its role in cosmology see |10j-|14j and 
references therein). The Casimir effect is common to all systems characterized by fluctuating 
quantities and has important implications on all scales, from cosmological to subnuclear. In 
the Kaluza-Klein-type models this effect has been used as a stabilization mechanism for moduli 
flelds which parametrize the size and the shape of the extra dimensions. The Casimir energy 
can also serve as a model for dark energy needed for the explanation of the present accelerated 
expansion of the universe (see [15] and references therein). In addition to its fundamental inter- 
est the Casimir effect also plays an important role in the fabrication and operation of nano- and 
micro-scale mechanical systems (see, for instance, jl6j ) and has become an increasingly popular 
topic in quantum field theory. 

The effects of the toroidal compactification of spatial dimensions on the properties of quan- 
tum vacuum for various spin flelds have been discussed by several authors (see, for instance, 
[l],[in]-[T^, [TT', T8] and references therein). In the present paper, we investigate one-loop 
quantum effects arising from vacuum fluctuations of a massive fermionic fleld on background 
of higher-dimensional spacetimes with an arbitrary number of toroidally compactified spatial 
dimensions. We will assume generalized periodicity conditions along compactified dimensions 
with arbitrary phases. Important quantities that characterize the quantum fluctuations are the 
fermionic condensate and the expectation value of the energy-momentum tensor. In the next 
section, by using the Abel-Plana summation formula, we derive a recurrence formula relating 
the fermionic condensates in topologies R^ x {S^^ and RP'^^ x (5^)*^"^. An alternative expres- 
sion for the topological part in the fermionic condensate is obtained by using the zeta function 
technique. In section [3] we consider the corresponding formulae for the vacuum expectation 
values of the energy-momentum tensor. In section [4] we give applications of general formulae 
to the Casimir effect for electrons in a carbon nanotube within the framework of 3-dimensional 
Dirac-like model. The main results of the paper are summarized in section [5l In appendix we 
show the equivalence of two representations for the vacuum expectation values obtained by the 
Abel-Plana summation formula and by the zeta function method. 

2 Fermionic condensate 

We consider a quantum fermionic field on background of {D + l)-dimensional flat spacetime 
with spatial topology x (S^)'^, p + q = D. The corresponding line element has the form 

D 



1=1 

where — oo < < oo, I = 1, . . . ,p, and ^ z'' ^ Li for I = p + 1, . . . , D. The dynamics of the 
field is governed by the Dirac equation 




(1) 



i'y^dfj^ip — mip = . 



(2) 
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In the (1) + l)-dimensional spacetime the Dirac matrices are N x N matrices with = 2[(^+^)/^] , 
where the square brackets mean the integer part of the enclosed expression. We will assume 
that these matrices are given in the chiral representation: 



with the relation ct^ct^ + cTpCT^ = 25^i,. For example, in = 4 the first four matrices 7^, 
^ = 0,1,2,3, can be taken the same as the corresponding matrices in 4-dimensional spacetime 
and 7"^ = 7''7^7^7^. In this case 0"i,cJ2,ct3 are the standard Pauli matrices and 

Note that, unlike to the Pauli matrices, (T4 is antihermitian. 

In this paper we are interested in the effects of non-trivial topology on the vacuum expec- 
tation values (VEVs) of the energy-momentum tensor and the fermionic condensate assuming 
that along the compactified dimensions the field obeys the boundary conditions (no summation 
over I = p + 1, . . . ^D) 

^{t, Zp,zg + Liei) = e^™' Zp, Zg), (5) 

with constant phases ai. In ([5]), Zp = (z^, . . . , z^) and Zg = {z^^^, . . . ,z^) denote the coordinates 
along uncompactified and compactified dimensions respectively, e; is the unit vector along the 
direction of the coordinate zK First we consider the fermionic condensate. 

For the topology under consideration the fermionic condensate (OlipiplO) (with |0) being the 
amplitude for the vacuum state) we will denote by {'ipip)p^q. We expand the field operator in 
terms of the complete set of positive and negative frequency eigenfunctions ''}: 

13 

where ap is the annihilation operator for particles, and 6^ is the creation operator for antipar- 
ticles. By using the commutation relations for these operators, the condensate is presented in 
the form of the mode-sum 

{i^i')p,, = Y.4~^(.x)i;'f\x). (7) 

/3 

In order to evaluate the condensate by this formula we need the explicit form of the eigenfunctions 
satisfying the boundary conditions ([5]). 

In accordance with the problem symmetry the dependence of these functions on the spacetime 
coordinates can be taken in the plane- wave form g*'' '""*'^*, uj = Vk"^ + rn?, with the wave vector 
k. From the Dirac equation we find 

(2P+i7rPy,a;)i/2 1^ (n • <t)u;W V^IT^ )' ^ ^ 



where (3 = (k, cr), n = k/fc, and a = (ui, (T2, . . . , 0"d), Vq = Lp+i-'-L/j is the volume of the 
compactified subspace. In these expressions a = 1, . . . , N/2, are one-column matrices 
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having Njl rows with the elements wf^ = 6ia-, and wi- ^ = iw^\ The eigenfunctions ([9]) 
are normahzed in accordance with the condition 

d^xV'f^+^^t^ = V- (10) 
In the discussion below we will decompose the wave vector into components along the uncom- 



pactified and compactified dimensions: k = (kp,kg), k = yk^ + k^. The eigenvalues for the 
components along the compactified dimensions are determined from the boundary conditions 

kg = (27r(np+i + Q;p+i)/Lp+i, . . . , 27r(nD + aD)/LD), rip+i, . . . , = 0, ±1, ±2, (11) 

For the components along the uncompactified dimensions one has —oo<ki < oo, I = 1, . . . ,p. 

Substituting the eigenfunctions Q into formula ([7]), for the fermionic condensate we find 
the expression 

with rig = (rip+i, . . . , n/j) and 

D 



g + ^ [27r(n, + aO/Li]2 + m^. (13) 
l=P+i 



We implicitly assume the presence of a cutoff function in (jl2p which makes the inegrosum finite. 

For the further evaluation of formula p2|) we apply to the sum over np+i the Abel-Plana 
summation formula in the form 1191 

(14) 



/-oo /-oo 

/(|np+i + ap+i|) = 2 / dxf{x)+i dx 



np+i=-oo ■'^ A=±l 



As a result, the fermionic condensate is presented in the decomposed form 

(V'^>p,<? = (V'V')p+l,g-l + Ap+i(Vi^/;)p,g, (15) 

where corresponds to the first term on the right-hand side of (|14p and is the fermionic 

condensate for the topology R^"*"^ x (S^)''"^. The second term on the right-hand side of formula 
(jlSp is induced by the compactness of the z^~^^ direction and is given by the formula 



II I 9 n,_iGZ9-i A=±l -"^"^-i 



2-1-^miV Lp+i ^ ^ p 

+i)/2r((p- 

where nq_i = (np+2, . . . , jt-d) and 



,2 



^ = [27r(n, + aO/L,]2 + m2. (17) 



Note that the expression on the right-hand side of (jl6l) is finite and the introduction of the cutoff 
function is necessary in the first term on the right of (jlSp only. 



4 



Expanding the function l/(e^ — 1) in the integrand of formula (jl6|) . we find an alternative 
form 



Ap+i(-^^)p,g = - (2^)p/2+iY ^cos{2Trnap+i) ^ ujl^_Jp/2{nLp+iuJn^_^). (18) 

V 1 n=l n,_ieZ9-i 

with the notation fu{x) = K^{x) / x'^ . From here it follows that in the case of periodic boundary 
condition along the direction z'^^^ (op+i = 0) the contribution to the fermionic condensate 
due to the compactness of the corresponding direction is always negative independently of the 
boundary conditions along the other directions. In the limit when the length of one of the 
compactified dimensions, say z\ I ^ p + 2, is large, the main contribution into the sum over ni 
in (lisp comes from large values of n/, and in the leading order we can replace the summation 
by the integration in accordance with 



1 +°° 1 t-c 

- V f{2^\ni + ai\/Li)^- / 



dyfiy)- 



n;=— oo 

The integral over y is evaluated by using the formula 



- / dy{y' + 62)^/Vp/2(cx/^^) = -=/(p+i)/2(c6), (19) 
Jo v27r 



and from (llSp the corresponding formula is obtained for the topology R^^^ x [S^)^^^. In the 
limit Li <^ -Z^p+i, I = p + 2, . . . , D, the main contribution into the topological part (jlSp comes 
from the term with nq_i = and in the leading order we have 



Ap+i(ViV)p,g — (2'Ky/'^+W^ ^ cos(27rnap+i)/p/2(nLp+im). (20) 



As we could expect, for large masses, mLpj^i ^ 1, the fermionic condensate given by formula 
p8|) is exponentially suppressed. 

After the recurring application of formula (|18p , the topological part of the fermionic conden- 
sate for spatial topology x (S^)^ is presented in the form 



D-l 



{i^i^)p,q = ^j + ii^^)3,D-j- (21) 



3=P 

For a massless field the fermionic condensate vanishes. 

An alternative form for the topological part in the fermionic condensate is obtained by 
making use of the zeta function technique [HI [2U]. We introduce the zeta function density 

^ ^ ' n,eZ'3 

with uj defined by relation p^ . In the case a/ = 0, m = 0, the point = is to be excluded 
from the sum. After the integration over kp, this function is presented in the form 



^ (47r)P/2r(s)F^ ^ 



D 



p/2-s 

i=P+i 
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An exponentially convergent expression for the analytic continuation of the function (|23|) is 
given by the generalized Chowla-Selberg formula [21]. The application of this formula to Eq. 
gives 

C(s) = Cm(s) + Cm(s), (24) 

where 

. ( f dkp 1 _ m^-^^ r{s~D/2) 

' J {2Tr)D (kl+m^Y (4^)^/2 T{s) ' ^ ' 

is the corresponding zeta function in the usual Minkowski spacetime and the part 

CpA^) = (2Tr)D/2Y( ) Yl cos{2TTuig-ag)fn/2-sim9(Lg,mg)), (26) 

with Lg = (ip+i, . . . ,L£)) and ctg = (op+i, ...,«£)), is induced by the nontrivial topology. The 
prime on the summation sign in ()26p means that the term nig = should be excluded from the 
sum and we have used the notation 




<7(L„m,)= >^ • (27) 



The topological part in (I24p is an analytic function at the physical point s = 1/2 and for the 
fermionic condensate one directly finds 

{tpip)p,q = ^Cp,g(l/2) = ~ (2t^){D+i)/2 2-. cos(27rmg • cXg)f(^D^iy2img(Lg,mg)). (28) 

In the case p = D — 1, q = 1 this formula coincides with (jlSp . In appendix we prove the equiv- 
alence of two representations (f2T]l and (f28l) for the topological part in the fermionic condensate 
for general case. Note that in (|28p we can write the function cos(27rmg • cxg) in the form of the 



product Hi^p+i cos(27rm.i 



3 Energy-momentum tensor 

In order to find the VEV for the operator of the energy-momentum tensor, we substitute the 
expansion ^ and the analog expansion for the operator tp into the corresponding expression 
for spinor fields, 

r^,{^, i;}='- [h(pd,)i' - (5(^^)7.) V-] ■ (29) 

Similar to the case of the fermionic condensate, by making use of the commutation relations for 
the annihilation and creation operators, one finds the following mode-sum formula 

(0|T^,|0) = {T^,)p,g = Y,T^Ai^'f\x),i;'f\x)} . (30) 

Substituting the eigenfunctions Q into this mode-sum formula, for the energy density and 
vacuum stresses one finds (no summation over I = 1, . . . , D) 

AT r 



6 



As in the case of the fermionic condensate, we will assume that some cutoff function is present, 
without writing it explicitly. 

After the application of summation formula (|14|) to the series over Up^i, we receive the 
following recurrence relation 

i'^Pp,q = (^/i')p+i,g-i + ^P+i{'^Pp,q^ (33) 

where (r^)p+i^g_i is the VEV of the energy-momentum tensor for the topology R^^^ x (S^y^^. 
The part Ap+i (T^)p_g is induced by the compactness of the z^^^ direction and is given by the 
expression (no summation over I) 

' - r((p + l)/2)y, ^2^^ L^_^ ,..,..-.2.^..,. _ 1 ' (34) 

with the notations 

n\ 4(n2 -Lo'i ^ 

= ^ ^12^, / = o,l,...,p, 

p+l 

f(p+^){u) = -2u\ f^'\u) =kf,l=p + 2,...D. (35) 
Expanding the integrand, this expression can also be presented in the form (no summation over 



INT °° 
\+iiTl)p,q = w2+'iV E E cos(2vrna,+i)<2^F(')(nL,+ia;„,_ J, (36) 



with the notations 



fW(^) = F^^\z)=fj,/,+,{z), l = l,...,p, 
piP+i)^,^ = -/p/2(^)-(p+l)W(^), (37) 
F«(z) = {ki/u;^^_,ffp/2{z), l=p + 2,...,D. 



It is easy to check that for a massless field the topological part (p6]) is traceless. As we see 
the vacuum stresses along the uncompactified dimensions are equal to the energy density. Of 
course, this property is a direct consequence of the boost invariance along the corresponding 
directions. In particular, from ([36]) it follows that in the case of periodic boundary conditions 
along the coordinate z^^^ (ctp+i = 0), the compactification along this coordinate increases the 
vacuum energy density independently of the boundary conditions along the other directions. The 
limiting cases of general formulae for the VEV of the energy-momentum tensor are investigated 
in a way similar to that described before for the condensate. 

From (j33p . for the VEV of the energy-momentum tensor in the topology x (S^)'^ one finds 

D-l 

{TPp,q=^^j+i{TPj,D-j. (38) 
j=p 

Now, by using the standard relations for the Mac-Donald function, it can be seen that the vacuum 
energy density and stresses along the compactified dimensions are related by the formula (no 
summation over /) 

dM{Tl^)p,,) = ^{Tl),^g,l=p + l,...,D. (39) 
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For the simplest Kaluza-Klein-type model with spatial topology x S^, from 
energy density one finds (ip+i = L, Op+i = a) 



for the 



3,1 



^ CO 

n=l 



cos{2ima) 



[(nmL)^ + 3nmL + 3]. 



(40) 



This quantity is positive for an untwisted field (a = 0) and is negative for a twisted field 
(a = 1/2). In the general case, the Casimir energy density is not a monotonic function of the 
size of the compactified dimension. This is seen from the left panel of figure [1] where we have 
plotted the quantity (|4U|) as a function of the parameter mL for different values of the phase a 
(numbers near the curves). The values of the phase are chosen in a way to show the transition 
from the positive energies to negative ones. In the right panel of figure [T] we have presented the 
Casimir energy density (|40p for a massless field as a function of the parameter a. 



V 





mL 



Figure 1: The Casimir energy density in the Kaluza-Klein-type model with spatial topology 
X 5"^ as a function of the parameter mL for different values of a (left panel). The right panel 
presents the corresponding quantity for a massless field as a function of a. 

An alternative expression for the VEV of the energy density is obtained by using the integral 
representation of the corresponding zeta function given by (j26p : 

m„eZ9 



{To)p,q = -YCp,g(-V2) = ^27r)(J'+i)/2 ^ cos(27rm, • Q!g)/(£,+i)/2(m5(L„m,)). (41) 



The equivalence of the representations (j38p and (j4ip for the energy density is seen in a way 
similar to that used in appendix for the fermionic condensate. The corresponding formulae 
for the vacuum stresses along compactified dimensions are obtained from relations ([39]) (no 
summation over I): 

P'l ~ (27r)(-D+i)/2 Yl "^?cos(27rmg •Q5)/(B+3)/2(m5(Lg,mg)), (42) 

with I = p + 1, . . . , D. A number of special cases of formula (j4ip for the Casimir energy can be 
found in literature (see [l].|10j-p^. |17j). For a massless fermionic field from (j4ip we find (no 
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summation over I) 

_ ,, r((Z) + l)/2) ^/ cos(27rm, • cx,) 

where I = p + 1, . . . ,D. Note that for a massless field the representation (|36|) has stronger 
convergence than the one given by (HSl), (01]): the summand in ([36j) decays exponentiahy instead 
of power-law decay in 



4 Applications to nanotubes 

In this section we specify the general results given above for the electrons on a carbon sheet 
rolled into a cylinder or torus making use of the description of the electronic states in terms 
of Dirac fermion fields. In this case D = 2 and we consider the geometries of cylindrical and 
toroidal nanotubes separately. Note that the Dirac-like model for electrons in a carbon nanotube 
is valid provided that the cylinder circumference is much larger than the interatomic spacing. 
For typical nanotubes the corresponding ratio can be between 10 and 20 and this approximation 
is adequate [5l [6] . 



4.1 Cylindrical nanotubes 

A single wall cylindrical nanotube is a graphene sheet rolled into a cylindrical shape. For this 
case we have spatial topology x with the compactified dimension of the length L. Note 
that the carbon nanotube is characterized by its chiral vector = riyuai + m^a2, with n^o, m^„ 
being integers, and L = |C/j| = a\/ n'^ + + nwiriw. In the expression for the chiral vector, 
ai and a2 are the basis vectors of the hexagonal lattice of graphene and a = |ai| = |a2| = 2. 46 A 
is the lattice constant. A zigzag nanotube corresponds to the special case C/j = (n^,0), and a 
armchair nanotube corresponds to the case Ch = (n^,n^). All other cases correspond to chiral 
nanotubes. The electron properties of carbon nanotubes can be either metallic or semiconductor 
like depending on the manner the cylinder is obtained from the graphene sheet. In the case 
riw — rriyj = 2>qw, Qw £ Z, the nanotube will be metallic and in the case — niw 7^ 3g«, the 
nanotube will be semiconductor with an energy gap inversely proportional to the diameter. In 
particular, the armchair nanotube is metallic and the (n^„,0) zigzag nanotube is metallic if and 
only if is an integer multiple of 3. 

In order to see the boundary conditions along the compactified dimension, we note that for 
the (n^,m^) nanotube the phase factor in the wavefunction is in the form Q^i^i+ii^-w-mw)/^]^ ^ 
mi € Z, where ip is the angular variable along the compact dimension. From here it follows that 
for metallic nanotubes we have periodic boundary conditions {ai = 0) and for semiconductor 
nanotubes, depending on the chiral vector, we have two classes of inequivalent boundary con- 
ditions corresponding to a/ = tt/3 (n^ — = 3qw + 2) and a/ = 27r/3 (n^ — = Sq^ + 1). 
In the expression for the Casimir densities the phases a; appear in the form cos(27rna/) and, 
hence, the Casimir energy density and stresses are the same for these two cases. 

Using the tight-binding approximation it can be seen that the electronic band structure close 
to the Dirac points shows a conical dispersion E(k) = f_p|k|, where k is the momentum measured 
relatively to the Dirac points and vp represents the Fermi velocity which plays the role of speed 
of light. The corresponding low-energy excitations can be described by a pair of two-component 
Weyl spinors, which are composed of the Bloch states residing on the two different sublattices 
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of the honeycomb lattice of the graphene sheet. The corresponding Fermi velocity is given by 
vp = 3ta/2 {vp ~ lO^cm/s in graphene), where t is the nearest neighbor hopping energy. Below, 
in specifying the formulae from previous section for the case D = 2, we consider a massive spinor 
field to keep the discussion general. The formulae for a massless case, appropriate for carbon 
nanotubes, will be given separately. 

In the case D = 2, the general formula for the fermionic condensate from section [2] takes the 
form {N = 2, p = l,q = l,Vq = L, Lp+i = L, a^+i = a) 



— Tfi 

= -Sc^imL), (45) 



where we have defined 



Sa(,x) = cos(27rna) 
=1 
1 



n 

n=l 



^ln[l-2e-^cos(2TO) + e-2^]. (46) 

In a similar way, for the VEV of the energy-momentum tensor from (|36|) we find (no sum- 
mation over /) 

1 „~nmL 

it}) 1 1 = — -S" cos{2TTna)G^^\nmL) (47) 

n=l 



with the notations 



G(0)(z) = GW(z) = 1 + z, G(2)(z) = -(2 + 2z + z^). (48) 
In particular, for the energy density we have 

{nKi = :^,S^^\rnL), (49) 

where the notation 

5(0) (x) = ycos(2™Q)e-"^^^^tZ^, (50) 

n=l 

is introduced. In figure [2] we have plotted the function Sa\x) for different values of a (numbers 
near the curves). In particular, the Casimir energy density is positive for armchair nanotubes 
(periodic boundary conditions). 
In the case m = we have 

(To°)i,i = (T/)i,i = -i(r|)i,i = (51) 

In particular, S^^\o) = 1.202, S^^^^iO) = -0.902, and 5j/Jj(0) = -0.534. Note that the corre- 
sponding fermionic condensate vanishes. In carbon nanotubes we have two sublattices and each 
of them gives the contribution to the Casimir densities given by (I5ip . So, for the Casimir energy 
density on a carbon nanotube with radius L one has 

(^o)&^ = ^5W(0), (52) 

where the standard units are restored. Hence, we see that the topological Casimir energy is 
positive for metallic nanotubes and is negative for semiconducting ones. 
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1 2 3 4 5 

X 



Figure 2: The function Sa (x) from (|50|) for different values of the parameter a (numbers near 
the curves). 



4.2 Toroidal nanotubes 

For the geometry of a toroidal nanotube we have the spatial topology (5"^)^ with p = and 
q = 2. In this case from the general formulae for the fermionic condensate we find 



m ^ SaAmLj) 2m ^ ^ , , , ^-m^miL'i+m^L'i 



VPW)o,2 = > — > > cos(27rmiai) cos(27rm2a2) 



(53) 

where the function Sq,(x) is defined by (06]). 

For the energy density ant the vacuum stresses the corresponding formulae have the form 
(no summation over /) 

, 0\ S^^{mLj) 2 ^ ^ cos(27rmiai)cos(27rm2a2) 1 + m-fi'(L2,m2) , , 

\Jo/o,2 - 2^ Ta \- - y, 7 77 XN 377 X — A54j 

,tt2 ^^i ^rnT^iml^i exp(mg(L2,m2)) 9\^2.V^2) 

rrv' V - 2 2 3 + 3x + 



(r/)o,2 = (To°)o,2 -—EE cos(27rm,-a,)^' 

j=l,2 mj=l 

^ E E cos{2TTmiai)cos{2TTm2a2)Lfmf U=mg(L2,m2). (55) 

mi=l 7712=1 

with I = 1,2 and (7(L2, m2) = a/ m'fLf + m2L2. Alternative expressions for the topological parts 
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are obtained from formulae (|36|) and (j38|) . For a massless field we find 



/rrO\ _ ^^J* (Q) , 2 ^ cos(27rmiai) cos(27rm2a2) 

^^0)0,2 - 2^ ^^ + -2^ 2^ KLf + miL2)3/2 ' ^^^^ 

j=l,2 J mi=lm2=l V i i z ^1 



3 ^ ,Lfmf 



■10 1 ^omf 

j=l,2r7ij=l J J 

6 2 2^°^(2'^"'-i'^i) cos(27rm2a2) 

-- 2^/^^"^/ (m?Lf + miLi)5/2 ' ^^^^ 

mi=lm2=l \ 1 i z z/ 

In particular, it is of interest to see the difference of the Casimir densities between the toroidal 
(with radii L\ and L2) and cylindrical (with radius L2) geometries of the carbon nanotube. For 
the condensate this difference is directly given by formula (jlSp and one has 

(V5V')o,2 = (^V')i,i-^Vcos(27rnai) V K^{n{UlL2)J^^\n2-.a2Y-.m^Ll^. (58) 

ra=l 712=— 00 

The first term on the right of this formula is the condensate for the topology i?^ x with the 
length of the compactified dimension L2 ■ Similar formula for the VEV of the energy- momentum 
tensor follows from (|36p (no summation over /): 

2 00 +00 
(T/)o,2 = (r/)i,i + --^ J^cos(2^nQi) J2 [4vr2(n2 + a2)^ + m2L^] 

2 n=l n2=— 00 



xF«(n(Li/L2)^47r2(n2 + a2)2 + m2L2), (59) 

where the functions F^''^ (z) are given by expressions (j37p with p = 0. The second terms on the 
right-hand sides of formulae (j58p and (|59p are induced by the compactification of the cylinder 
(with radius L2) along its axis. In figure [3] we have plotted these terms for the energy density, 
Ai(Tq)o,2 (left panel), and for the stress along the axis of the cylinder, Ai(Tj^)o,2 (right panel), 
for a massless fermionic field as functions of the ratio L1/L2. The numbers near the curves 
correspond to the values of (01,02)- As we have mentioned before the values of the phase 
ai = 0,1/3 are realized in carbon nanotubes. The vacuum stress Ai(r2)o,2 is related to the 
quantities plotted in figure [3] by the zero trace condition for the energy-momentum tensor of a 
massless field. 

The corresponding formulae for the Casimir densities in toroidal nanotubes, which we denote 
by {TI)q^ , are obtained from (j56p . (|57p and (j59p in the massless limit with additional factor 2 
which takes into account the presence of two sublattices: (T/)J2^ = 2(T/) 

o,2|m=o- In standard 

units the factor hvp appears as well. Note that if the chiral vector Ch is directed along the axis 
z"^ then one has L2 = \Ch\- The translational vector defining the unit cell, T, is perpendicular 
to Ch and its components are related to the components of the chiral vector by the formula 

^ = dr d^^^' ^^'^ 

where dr = gcd{nyj , rriyj) if (m^ — n^) is not a multiple of 3 x gcd(n^t,, m^^,) and = 3 x 
gcd{nw,mw) if (m^ — n^) is a multiple of 3 x gcd{nw,mw)- Here gcd means the greatest 
common divisor. Now for the length of the second toroidal dimension we have Li = A^^|T|, 
where is the number of unit cells along the corresponding direction. By taking into account 
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Figure 3: The difference between the vacuum energy densities (left panel) and stresses (right 
panel) between the cylindrical (with radius L2) and toroidal (with radii Li and L2) geometries 
for a massless fermionic field. The numbers near the curves are the corresponding values for 
(ai, 02). 

that |T| = \/3Li/(ir, for the ratio of the lengths of the torus in ()59p one finds L1/L2 = V^Nw/dr. 
From the graphs in figure[3]it follows that the toroidal compactification of a cylindrical nanotube 
along its axis increases the Casimir energy for periodic boundary conditions (ai = 0) and 
decreases the Casimir energy for the semiconducting-type compactifications. In particular, the 
Casimir energy of the armchair cylindrical nanotube increases by the compactification if is 
an integer multiple of 3 and decreases otherwise. 

5 Conclusion 

In the present paper we have investigated the topological Casimir effect for a massive spinor 
field on background of spacetime with an arbitrary number toroidally compactified spatial di- 
mensions. The boundary conditions along compactified dimensions are taken in general form 
with arbitrary phases. For the evaluation of the Casimir densities we have used the direct 
mode-summation method. By applying to the corresponding mode-sums the Abel-Plana for- 
mula, we have derived recurrence formulae which relate the VEVs for the topologies x {S^ Y 
and BP^^ x {S^Y~^ . The part induced by the compactness of the {p + l)-th direction is given 
by expression (jlSp for the fermionic condensate and by expression (|36p for the VEV of the 
energy-momentum tensor. The total topological VEVs are obtained after the summation over 
all compactified dimensions, formulae (I2ip and (I38p . Alternative expressions are obtained by 
using the generalized Chowla-Selberg formula for the analytic continuation of the corresponding 
zeta function. These expressions are given by formula ()28p for the condensate and by formulae 
(|4ip and (j42p for the energy density and vacuum stresses along compactified dimensions. Note 
that the stresses along the uncompactified dimensions coincide with the energy density. This 
property is a direct consequence of the boost invariance along the corresponding directions. For 
a massless fermionic field the condensate vanishes and the expressions for the VEVs of the en- 
ergy density and vacuum stresses take the form ()43p and (j44p . Note that, unlike to the case 
of a massive field, the convergence of the multiseries in the latter case is power-law. In the 
representation based on the application of the Abel-Plana summation formula we have expo- 
nentially convergent multiseries in both cases of massive and massless fields. On the example 
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of the simplest Kaluza-Klein-type model with spatial topology x we have demonstrated 
that, unlike to the special cases of twisted and untwisted fields, in general, the Casimir energy 
density is not a monotonic function of the size of the internal space. 

In section m we specify the general formulae for the model with D = 2. This model may be 
used for the evaluation of the Casimir densities within the framework of the Dirac-like theory 
for the description of the electronic states in carbon nanotubes where the role of speed of light 
is played by the Fermi velocity. Though the corresponding spinor field is massless, to keep the 
discussion general, we present the formulae for the cylindrical and toroidal geometries in the 
massive case and specify the results for the nanotubes separately. For carbon nanotubes the 
fermionic condensate vanishes and the VEV of the energy-momentum tensor is given by formulae 
(I52p for cylindrical nanotubes and by ()56p and (j57p (with an additional factor 2 which takes into 
account the presence of two sublattices) for toroidal nanotubes. In the case of toroidal nanotubes 
an alternative representation with the stronger convergence of the series is given by formula 
(I59p with 171 = 0. The topological Casimir energy is positive for metallic cylindrical nanotubes 
and is negative for semiconducting ones. We have shown that the toroidal compactification of a 
cylindrical nanotube along its axis increases the Casimir energy for periodic boundary conditions 
and decreases the Casimir energy for the semiconducting-type compactifications. In particular, 
the Casimir energy of the armchair cylindrical nanotube increases by the compactification if the 
number of unit cells along the axis of cylinder is an integer multiple of 3 and decreases otherwise. 
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A Equivalence of two approaches 

In this section we show that the formulae (j2ip and (128|) for the topological part in the fermionic 
condensate are equivalent. First of all we note that from formula (j28p one has 




X 





(61) 



mq_ieZ9-l 



Hence, we should prove the relation 




(62) 



For this we will use the Poisson's resummation formula 




(63) 



for the function 



F{x) = cos(27rx • Q;g_i)/(£,_i)/2(mW c/2(l^_^^ x) -h 



'p+i"^p+i)- 



(64) 
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After the integration over x we find 

^ COs(27rm,_i • OLq^i)f(^D-l)/2{mg{'Lq, Ulq)) 



= ^ / dx cospTTX • {(Xq^i + ng_i)] /(£,_i)/2 (m y^c/^ (Ly_i , x) + Ll^^ml^^).{65) 

ng_ieZ9-i 

For the evaluation of the integral on the right hand side we first introduce a new integration 
variables in accordance with yi = XiLi and then introduce spherical coordinates. The integration 
over the angular coordinates is expressed in terms of the Bessel function. At the final step the 
integral is evaluated by using the formula |22j 

dyy^^^J,{by) /.(c^/^T^) = ^{b" + c^Y'^-^ U-,-i{aV^^). (66) 
This leads to the following result 



j dx cos(27rx • (ag^i + r^q^l))f(o-l)/2{'m^Jg'^{"Lq^l,yi) + L. 



(27r)('^-^)/^L,+i „ 



m^-Wq "'^n9-i-^P/2("^P+i^P+i^n,_i), (67) 

where Wn^.i is defined by relation (fT3]) . Substituting this relation into ([65]) leads to the result 
which proves the equivalence of two expressions for the topological part. 
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